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A general class of wireless relay networks with a single source-destination pair is considered. 
Intermediate nodes in the network employ an amplify-and-forward scheme to relay their input signals. 
In this case the overall input-output channel from the source via the relays to the destination effectively 
behaves as an intersymbol interference channel with colored noise. Unlike previous work we formulate 
the problem of the maximum achievable rate in this setting as an optimization problem with no 
assumption on the network size, topology, and received signal-to-noise ratio. Previous work considered 
only scenarios wherein relays use all their power to amplify their received signals. We demonstrate that 
this may not always maximize the maximal achievable rate in amplify-and-forward relay networks. The 
proposed formulation allows us to not only recover known results on the performance of the amplify- 
and-forward schemes for some simple relay networks but also characterize the performance of more 
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complex amplify-and-forward relay networks which cannot be addressed in a straightforward manner 
I> using existing approaches. 

Using cut-set arguments, we derive simple upper bounds on the capacity of general wireless relay 
networks. Through various examples, we show that a large class of ampUfy-and-forward relay networks 
can achieve rates within a constant factor of these upper bounds asymptotically in network parameters. 



I. Introduction 

Since their introduction in [[l| Amplify-and-Forward (AF) relay schemes have been studied in 
the context of cooperative communication [|2j, [[3j, estimating the capacity of relay networks [4|- 
[|6|, and analog network coding [|7|-pT|. For cooperative communication, AF schemes provide 
spatial diversity to fight against fading; for capacity estimation of relay networks, such schemes 
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provide achievable lower bounds that are known to be optimal in some communication scenarios; 
and for analog network coding, given the broadcast nature of the wireless medium that allows the 
mixing of the signals in the air, these schemes provide a communication strategy that achieves 
high throughput with low computational complexity at internal nodes. In this paper, we concern 
ourselves mostly with the capacity analysis of a general class of Gaussian AF relay networks. 
Extensions of our method and results to cooperative communication and analog network coding 
scenarios is part of our future work. 

In previous work, while analyzing the performance of AF schemes in relay networks one or 
more of the following assumptions have been made: networks with a small number of nodes 
|(8|, pO| ; networks with simple topologies [|3|-[|5|, [[8|, pO| ; or relay operation in the high- 



SNR regime, [10|. However, for two reasons, we believe that it is important to characterize the 
performance of the AF schemes without such assumptions. First, we feel that for a scheme such 
as amplify-and-forward that allows one to exploit the broadcast nature of the wireless medium 
such assumptions on the size and topology may result in lower achievable performance than 
otherwise. Second, even in the low-SNR regimes amplify-and-forward can be capacity-achieving 
relay strategy in some scenarios, [5]. Therefore, a framework to address the performance of AF 
schemes in general wireless relay networks is desired. 

However, one major issue with constructing such a framework is the following. In general 
wireless relay networks with AF relaying, the resulting input-output channel between the source 
and the destination is an intersymbol interference (ISI) channel (Q, pO| ) with colored noise. 
This is because both the source signal and the noise introduced at the relay nodes may reach the 
destination via multiple paths with differing delays. Without the assumptions above, this results 
in a formidable problem to analyze with the existing methods [|4|. 

Our main contribution is that we provide a framework to compute the maximum achievable 
rate with AF schemes for a class of general wireless relay networks, namely Gaussian relay 
networks. This framework casts the problem of computing the maximum rate achievable with 
AF relay networks as an optimization problem. Our work shows that amplifying the received 
signal to the maximum possible value at intermediate nodes might result in sub-optimal end-to- 
end throughput. Also, we establish the generality of the proposed formulation by showing that 
it allows us to derive in a unified and simple manner not only the various existing results on 
the performance of simple AF relay networks but also new results for more complex networks 
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that cannot be addressed in a straightforward manner with existing methods. We show through 
various examples that for a large class of relay networks the AF schemes can achieve rates 
within a constant factor of the cut-set upper-bounds on the capacity of general wireless relay 
networks. 

The paper is organized as follows. In Section |ll] we introduce the general class of Gaussian AF 
relay networks addressed in this paper. In Section Ini] we formulate the problem of maximum rate 



achievable via AF schemes in these networks. In Section IV we compute the rates achievable via 
AF schemes for two instances of such relay networks under various communication scenarios, 
and then in Section |V] we discuss the asymptotic behavior of the gap between these rates and 
the corresponding upper bounds on the capacity of general wireless relay networks computed 



there. Section VI concludes the paper with a summary. 



II. System Model 

Let us consider a (M + 2)-node wireless relay network with source s destination t and M 
relays as a directed graph G = (V, E) with bidirectional links, as shown in Figure [T] Each node 
in the network is assumed to have a single antenna. Let us assume that the degree of the source 
node is A^^ + 1, with it being connected to the destination node and a subset Sg of the relay 
nodes, A^^ = l^sl. Similarly, let us assume that the degree of the destination node is Nt + 1, with 
it being connected to the source node and a subset St of the relay nodes, Nt = \Ss\. In general, 
SsUStCV\{s,t}. 



4 



At instant n, the channel output at node i,i &V \ {s}, h 

ViVA = 5^ hjiXj[n\ + Zi[n\, — oo < n < oo^ (1) 

where Xj[n] is the channel input of the node j in the neighbor set A/'(i) of node i. In ([T|), 
hji is a real number representing the channel gain along the link from relay j to relay i. It is 
assumed to be fixed (for example, as in a single realization of a fading process) and known 
throughout the network. Further, is a sequence (in n) of independently and identically 

distributed (i.i.d.) Gaussian random variables with zero mean and variance a^, Zi[n] ~ A/'(0, a^). 
We also assume that Zi are independent of the input signal and of each other. The source symbols 
Xs[n], — oo < n < oo, are i.i.d. Gaussian random variables with zero mean and variance Pg that 
satisfies an average source power constraint, Xs[n] ~ J\f(0,Ps). We assume that the i"' relay 
node's transmit power is constrained as: 

Efx- [n]] < Pi, -oo < n < oo (2) 

In a general wireless relay network there may exists cycles. If a relay merely amplifies and 
forwards its received signal in such scenarios then it may be use a significant fraction of its power 
budget on forwarding the previously forwarded information. Therefore, motivated by some work 
on analog network coding, such as [12], we propose the following relay operation to allow the 
relays to expend their transmit power in forwarding only the "new" information. 

Relay operation: We assume that each relay node maintains a buffer of signals it forwarded 
previously. Therefore, each relay node i, after receiving the channel output yi[n] at time instant 
n executes the following series of steps: 

Step 1: Obtain the residual signal ?/-[n] from its input yi[n] by subtracting the contributions 
of previously forwarded signals (if any) from yi[n]. 

Step 2: Compute the power ^ of the residual signal y[ [n] . 

Step 3: At instant n + 1 transmit the scaled version of the residual signal y • [n] of its input at 
time instant n: 

x,[n + l] = (3,y'M, ^<Pl<Plma. = P^/P'R,^ (3) 

where the I3i is the scaling factoiQ Let the network-wide amplification vector for the M relay 

'Note that, in general, Pi may depend on relay's past observations. I3i[n\ — }i,n{Yi[n — 1], . . . ,y4l]). However, due to 
practical considerations, such as low-complexity operation, we do not consider such scenarios here. 
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nodes be denoted as f3 = . . . , (3m)- 

Remark 1: One of the major advantages of this relay operation is that by subtracting the 
previously forwarded signal from its input, a relay node expends its power in forwarding only 
the "new" information. In a general wireless relay network there may exists cycles. In such 
scenarios if a relay merely forwards its received signal then it may result in the relay forwarding 
the scaled version of the linear combination of two or more of its previously forwarded signals. 



This may provide higher achievable rate as in [13| 



Using ([T]) and ([3]), the input-output channel between the source and destination can be written 
as an intersymbol interference (ISI) channel that at instant n is given by 



yt[n] = hstXs[n] + zt[n] 



d=l {ii,...,id)€Kd 



Xs [n — d] 



zi [n — d] 



(4) 



+ ^ (^AlhMh ■ ■ ■ hid-iidf^idhidt 

d=l {ii,...,id)&KM.d 



ZM[n - d], 



where Kd-, 1 < d < D^, is the set of (i-tuples of node indices corresponding to all paths from 
the source to the destination with path delay d and is the length of the longest such path. 
Note that along such paths < M . Similarly, Km,d, 1 < m < M,l < d < D"^, is the set of 
d-tuples of node indices corresponding to all paths from the relay to the destination with 
path delay d, is the length of the longest such path from m"^ relay to the destination. It 
should be noted that max{D^, . . . , D^^) = - 1. 

Let us introduce modified channel gains as follows. For all the paths between the source s 
and the destination t: 



ho — hgt 

hd = hsi^fii-^^hi-^i^ . . . hi^_-^i^l3i^hi^f, ^ ^ d < 

{il,...,id)&Kd 



(5) 



1 2 



M-1 M 



Fig. 2. Type A Relay Network. 

For all the paths between the m*-relay, 1 < m < M, and the destination t: 

hm,o = 0, (6) 

hm,d ^ ^ Pmhmii ■ ■ • ^*d-i*d/^*d^*dt' ^ d D 

Remark 2: Note that though there may be exponentially large number of paths between the 
source and the destination as well as between a relay and the destination, the modified channel 
gains for all such paths as in ([5]) and ([6]) can be efficiently computed using the line-graphs [ |14J . 

In terms of these modified channel parameters, the source-destination ISI channel in Q can 
be written as: 

yt l'^] = X] ^i^^ ~ -^l + 5Z N - i] + • • • + hujZM [n - j] + zt [n] (7) 

i=o i=o j=o 

Before proceeding further, let us introduce two special cases of the general class of relay 
networks introduced in the beginning of this section. In the rest of this paper, we illustrate 
various concepts and derive some results using these two special networks. 

Type A Network: For one source-destination pair and M relays. Type A network is defined 
as: G = {V,E), where V = {s,t,l, . . . , M} and E = {{s,t), {s,i), {i,t) : i G {1,...,M}}. 
In other words, the source node shares an edge with the destination node and each of M relay 
nodes. Similarly, the destination node shares an edge with the source node and each of M relay 
nodes. However, no pair of relay nodes share an edge between themselves. Figure |2] illustrates 
an instance of Type A networks. 

Type B Network: For one source-destination pair and M relays. Type B networks are defined 
as: G = iV,E), where V = {s, t, 1, . . . , M} and E = {{s,t), {s,i), {i,t), {j, j + 1) : i e 



Fig. 3. Type B Relay Network. 



{1, . . . , M},j G {1, . . . , M — 1}}. In Other words, as in Type-A networks, the source node 
shares an edge with the destination node and each of M relay nodes. Similarly, the destination 
node shares an edge with the source node and each of M relay nodes. However, relay shares 
an edge with j + node, for all j G {1, . . . , M — 1} . Figure [3] illustrates an instance of Type 
B networks. 



III. Achievable rates for the source-destination ISI channel in general AF 

Relay Networks 

We first derive the expression for the achievable rate for the source-destination channel in Q 
for a given amplification- vector /3 and then formulate the problem of maximizing the achievable 
rate over the domain of feasible /3. 

Lemma 1: For given length-M vector /3, the achievable rate for the channel in (|7]) with i.i.d. 
Gaussian input is: 



1 

2^ 



log 



1 + 



where 



H{X) = J2hje-'^\ Hm{X) = Y.hm,,e- 



ijX 



(8) 



(9) 



j=0 



j=0 



Proof: In [ 15 1 a Discrete Fourier Transform (DFT) based formalism is developed to compute 
the capacity of Gaussian channel with source intersymbol interference (ISI). We compute the 
maximum achievable rate for the channel in (|7]) for a given amplification- vector f3 by generalizing 
this formalism to also include the ISI channel for the Gaussian noise at each relay node resulting 
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in colored Gaussian noise at the destinatiorj^ The details of the proof are in the Appendix [a] 

■ 

Remark 3: The derivation of an expression for I{Ps, l3) with jointly Gaussian inputs is similar 
to the proof of Lemma [T| However as such an expression does not aid in the presentation of 
our ideas we do not discuss it in this paper. 

For a given network-wide amplification vector /3, the achievable information rate is given by 
I{Ps,/3). Therefore the maximum information-rate Iaf{Ps) achievable in an AF relay network 
with i.i.d. Gaussian input is defined as the maximum of I{Ps,i3) over all feasible f3, subject to 
per relay-node amplification constraint ([3]). In other words: 

(PI): Iaf{Ps)= max J(P„/3) (10) 



Substituting for H{\) and Hm{\) from ^ in ([8]), we can rewrite problem ( fTO] ) equivalently as: 



(P2):/af(P.) = max /(Ps,/3) (H) 



nPs,f3) = -L / log 







dX 



^'Ef=>^cos(zA). 

The coefficients At, 1 < i < D^, and Bi, 1 < i < -D", are defined in Appendix |B] 

The formulation of the problem P2 is illustrated in Appendix |C] for the Type A and Type B 
relay networks introduced in the previous section. 

In general AF relay networks the simultaneous relay transmissions may interfere and if the 
relays always amplify the received signals to the maximum possible then it may result in sub- 
optimal end-to-end throughput. Therefore the scaling factor for each relay must be optimally 
chosen to maximize the achievable rate. This is emphasized by PI and its significance is 
illustrated by the following example. 

2 n n 

Note that the capacity of discrete-time Gaussian ISI channel was known prior to |15|, for example in 1251-1271. However, 
the analysis in these papers, based on the asymptotic properties of Toeplitz forms (28| , (29| , is not easily amenable to derive 
the capacity results of the ISI channel with colored Gaussian noise. Therefore, like |30| , (31| , we too use more accessible DFT 
based formalism developed in f 15J to compute the maximum achievable rate for the ISI channel with colored Gaussian noise 
in 0. 
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Fig. 4. Diamond network. 





/?2 



Fig. 5. Tlie achievable rate for tlie Example 1 when /3i = /3i, max and /32 lies in [— /32,maa;, /32 



Example 1: Let us consider the relay network in Figure |4j Let hgi = l,hs2 = 0.1, /iit = 
/i2t = 1. Let Pg = Pi = P2 = 10 and noise variance cr^ = 0.1 at each node. Therefore, we have 

Pi 



hliPs + cy^ 
2 _ -P2 

2,max - ^2^p^ ^ ^2 



0.99 
50.0 
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In this case from (fTO]), we have the following rate maximization problem: 



Iaf = max - log 

/3i,/32 2 



1 + 100 



+ 0.1/32)' 



subject to constraints < Pf < (3l^^^ and < /3| < Pl^ax- 

The optimal solution of this problem is = 0.995, /32 = 0.225). The objective function is 
plotted in the Figure |5] for /3i = 0.995. Therefore, it follows that in this case (32 = (32,max is not 
the optimal amplification factor. ■ 

With this observation and the definition of the relay operation given in the previous section, it is 
appropriate to call the forwarding scheme proposed in this paper as subtract-scale-and-forward. 

Remark 4 (the computational complexity of the problem P2): Writing the objective function of 



the problem P2 as sum of ratios, as in ( |A.14[ ), we can show from [16| that each of these ratios 
is neither quasiconcavej^ (therefore, more than one local optimum may exist) nor quasiconvex 
(therefore, the globally optimal solution may not exist at an extreme point or on the boundary of 
the domain of optimization). The sum of such non-quasiconcave and non-quasiconvex ratios is 
also non-quasiconcave and non-quasiconvex. Therefore, the problem P2 is a global optimization 
problem, specifically it belongs to a subset of global optimization problems, called sum-of-ratios 
problems, which are supposed be hard, in general, |17|. Further, in general, Pi^max, ^ <i < M, 

def 

depends on /3_j = . . . , /3j+i, . . . , (3m)- This dependence of (3i^max on f3^i makes the 
constraint set in ( [TT] ) non-convex for all but Type A relay networks. Based on these arguments, 
we conjecture the problem P2 to be computationally hard. However, at present we do not have 
a formal computational complexity-theoretic proof and constructing such a proof is the part of 
our future work. 

Remark 5 (the approximation schemes to solve the problem P2): The lack of the exact 



complexity classification of the optimization problem in ( [TT] ) notwithstanding, we envision 
solving this problem to be an important step towards characterizing the maximum rates achievable 
with various other relaying schemes in general wireless relay networks as such schemes also 



result in problem formulations similar to ( [TT] ). In [32|, we concern ourselves with developing a 
unified framework to efficiently approximate the optimization problems such as ( [TT] ). 



'a function / : 
then f{Xx + (1 - 



S ^ R defined on a convex subset S of real vector space is quasiconcave if whenever x,y £ S and A £ [0, 1], 
A)y)>min(/(x),/(y)), (l6j. 
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In the following we show that the problem formulation above not only allows us to reproduce 
existing results on the achievable rates for some special classes of amplify-and-forward wireless 
relay networks in a simple manner, but also allows us to compute the achievable rates for 
much broader class of amplify-and-forward wireless relay networks which could not be hitherto 
addressed with existing methods. 

IV. Approximating the Maximum Achievable Rate Iaf{Ps) for 
Amplify- AND-FoRWARD Relay Networks 

Let us consider the problem P2 when the relays operate instantaneously as in [8], that is, the 
relays amplify-and-forward their input signals without delay. Therefore, for relay node i,l < i < 
M, we have Xi[n] = (3iy'i\n]. Note the possible system instability resulting from this assumption 
is avoided by the relay-operation (buffering and subtracting of the previously forwarded signals) 
as given in the Section |llj With this assumption, P2 reduces to (after setting A = in ( [TT] ) and 
then integrating): 



(93): lAFiPs) = max -log 

/3:0</32</32 2 



1 + 



(12) 



The problem P3 can also be derived directly from the channel model in (|7]), however we do not 
discuss that derivation here for the sake of brevity. 

Given the particular form of functional dependence of AjS and i?jS on the components of the 
amplification vector (3 (as follows from Q, ([5]), ( |B.2| ), and ( |B.9[ )) and dependence of /^jS on each 



other, the problem P3 is a Geometric Program (GP) [18|, where both, the objective function as 
well as the constraints are expressed in terms of posynomials. However, minimizing or upper 
bounding (equivalently, maximizing or lower bounding) a ratio of two posynomials is a non- 



convex problem that is intrinsically intractable [ 19 1, |20rl Therefore, for a general relay network. 



it is not possible to exactly solve or lower bound the problem P3 in a computationally efficient 



manner. However, unlike the problem P2 in ( [TT] ), for the problem P3, efficient approximation 



schemes exist that solve such problems iteratively by solving a series of GPs, as discussed in [ 19 



''Note that in (19| , (20| the objective function can be written as an inverted posynomial in high-SINR regime. Therefore the 
problem is efficiently solvable using GP methods. It is only in low to moderate SINR regimes where the objective function 
cannot be so written, the problem is NP-hard. However, the objective function in \\2\ involves the ratio of posynomials for all 
SNR values. Therefore, in general, the problem is hard irrespective of SNR value. 
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Section 3.3]. Also, for some specific relay networks, under some assumptions, we can efficiently 
compute the lower bounds, as we show next. 

Let us consider the problem P3 when the relay nodes are constrained to use the same 
amplification-factor, that is, (3i = (3, for all 1 < i < M. In the practical setting, this assumption 
considerably simplifies the system-design with (3 set to one particular value for all relay nodes. 
Then, (P3) reduces to 

(P4): lAFiPs)= max (13) 

^ — f-' —h'max 



lAF{Ps.P) = -\og 



1' I + 



(14) 



Note that the solution of (P3) cannot be smaller than the solution of (P4) because the set of 
feasible (3 : [5i = (5, for (P4) is a subset of the set of feasible (3 for (P3). 

Note 1: In general, ~ ^i- Therefore ^ (^fmax ~ S -^j- However, as we are considering 

the case of equal (3, jSi^rnax = Pmax, so we have Pi or ~ M'^ Yl Pi- 

Now we discuss solving the problem P4 for Type-A and Type-B relay networks introduced in 
Section |n[ in different communication scenarios. The proofs of various propositions and lemmas 
are provided in the Appendix |Dj 

A. Type-A Relay Network 

Let us first consider Type-A relay network as in Figure |2j For such networks, we solve the 
problem P4 in the following two scenarios. 

Scenario 1 (No attenuation network): Let us assume that there is no attenuation along any 
Unk in the network, that is, hgt = hsi = = 1 for all 1 < i < M. The problem P4 in this case 
is: 

Proposition 1: 



Iaf{Ps)= max J log 



-^^Ps{l + M(5f 



a2 1 + M^2 

Lemma 2: Iaf{Ps) attains its global maximum at (3opt = 1- 
Now let us consider two particular ways in which (3max varies with network size. 
Scenario 1, Case A (Increasing relay power): For the given network with M relay nodes, let 
us assume that the sum power of the relay nodes is constrained as follows: 

M AI 

E[X^] < ^ Pi < M"+iQ, u>0,Q = constant 

m=l m=l 
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So = M^Q. From Lemma |2| we have for M — )• oo. 

ilog[l + ^(l + M)], if/3™a. >1, 
Iaf{Ps) = < (15) 

I log[l + 4(M + -T^)], otherwise 

V 

Scenario 1, Case B ( Constant total relay power): Let us consider the case where the sum power 
of relay nodes is fixed irrespective of the number of relay nodes, that is J2m=i Pi ^ QiQ = 
constant. Therefore, we set = As M — )■ oo, for sufficiently large M, Pmax < Popt = 1- 
Therefore, from Lemma [2[ /3 = (3max maximizes the achievable rate and we have for Af — )■ oo 

1 



Iaf{Ps 



log 



^2 1 + g 



(16) 



Scenario 2 (Bounded channel gains): In Scenario 1, we considered no attenuation relay 
networks. Now, let us consider the scenario where the channel gains are arbitrary, but strictly 
bounded, < hst, hsi, hu < oo, 1 < i < M. The problem P4 in this case is: 

Proposition 2: 



Iaf{Ps) 



max - log 



Ps {hst + P Y^fii hsjh 



Lemma 3: Iaf{Ps) attains its global maximum at (3opt = T^'^kP .2 • 

Increasing relay power: Let us consider the increasing total relay power scenario as in Scenario 
1, Case A. Let = M^Q. In this case, following Lemma [s} we obtain the following lower 
bound on the achievable rate as M — )■ oo: 



Iaf{Ps) > 2 log 



if /3. 



Iaf{Ps) > ^log 



P. 



{Mhmin + 1 



h ■ 

^ "'s,max""rmn 
hm.n.x 



(17) 



> 



E£i hs.hu 



hst Ei=i hi ■ 



l+^{Mhr, 



+ 



h ■ 

, "'s,max'^min 



fin 



where h^ 



otherwise, 

max{hst, hsi, • • • , hsn}, h^in = min{/iit, . . . , liMt}, and hmax = max{hit, hut}- 
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B. Type-B Relay Network 

Let us consider the Type-B relay network as in Figure [3| For such networks, we consider the 
no-attenuation scenario where all channel gains are set to unity, that is, hgt = hgi = hu = 1, 
I < i < M as well as /ij,t+i = 1, 1 < i < M — 1. The problem P4 in this case is: 

Proposition 3: 

1 



Iaf{Ps,/3) 



log 
log 



1 + 



a 



1 + 



1 + E£i/3^(^^^^i^)' 



n2 



1 + 



[-3 



2^ 



'1 



^-1 



2/3^ + M] 



Iaf{Ps 



max 



(18) 



It can be proved that the objective function is quasiconcave, therefore a unique global max- 
imum exists. Let (3 that solves ( [T8| ) be denoted as /3 = /3op(. However, obtaining a closed- form 
expression for (3 opt does not appear straightforward, though it can be numerically computed for 
any M. 

Constant total relay power: Let us consider the case where the sum power of relay nodes is 
fixed irrespective of the number of relay nodes as in Scenario 1, Case B. Let P'^ax ~ m- 
M — )■ oo, for sufficiently large M, (3max < f^opt- Therefore, (3 = (3max maximizes the achievable 
rate and we have the following rate achievable asymptotically as M — )■ oo 

1 



lAFiPs 



log 



1 + 



QM 1 + ^Q/M' 



(19) 



V. Asymptotic Capacity 



In the previous two sections, we formulated the problem of maximum achievable rate for 
AF relay networks and then we computed explicit lower bounds to the capacity of two specific 
AF relay networks in various communication scenarios. In this section, we first derive an upper 
bound to the capacity of the general relay networks we address in this paper as introduced in 
the Section [Ilj We then discuss the asymptotic behavior of the gap between this upper bound 
and the lower bounds computed for two specific relay networks in the previous section. 
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Fig. 6. Illustration of the particular broadcast and multiple-access cuts of the Gaussian relay network in Figure [T] 

A. Upper bounds to Capacity of Relay Networks 

In Q, an upper bound to the capacity of the relay network of Type A is computed using 



a weaker corollary of the cut-set theorem [22 Theorem 15.10.1] and the capacity formula for 



Gaussian vector channels with fixed transfer function p3| . Using this corollary and the capacity 
formula, we can also compute the upper bound to the capacity of the general relay network 
introduced in the Section as stated in the following proposition. 

Proposition 4: The capacity C of a general wireless relay network is upper-bounded as C < 
minjCsc, Cmac}, where Cbc and Cmac are the upper bounds on the capacity of the broadcast 
cut and multiple-access cut respectively, as in the Figure |6} and are given as follows 

Ps 



C 



BC 



log 



C 



MAC 



log 



1 + 



Proof: The details of the proof are provided in the Appendix |Ej ■ 
Remark 6: Proposition 3 in [|4J can be obtained as a special case of the above proposition by 
setting Ns = Nt = M. 



B. Type A Relay Networks 

No attenuation, increasing relay power: In this case, the broadcast bound Cbc is always 
asymptotically smaller than the multiple-access bound Cmac^ as follows from 



^^mac—Cbc ~ X -|- 



-, for large M. 
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Therefore it suffices to compute the asymptotic gap between Cbc and the lower bound in ( [T5] ). 
In fact, in this case we have 

Cbc - 2Iaf{Ps) = 0, for all M > 1 

The actual capacity C of the relay network in this case is bounded by Csc/'^ < C < Cbc- 

No attenuation, constant total relay power: In this case also the broadcast bound Cbc is 
asymptotically smaller than the multiple access-bound Cmac, as shown below 

lim 2'^*^'^c'~^s^ = 1 + ^ 

M^oo P 



Therefore we only address the asymptotic gap between Cbc and the lower bound in ( [T6| ). We 
have 

lim Cbc - 2/af(Ps) = J log(l + l/Q) 

The actual capacity C of the relay network in this case is bounded by | (Cbc* — ^ log(l + l/Q)) < 
C < Cbc and the bound gets tighter with increasing Q. 



Bounded channel gains: The gap between Cbc and lower bound of achievable rate in (IT) is 
bounded asymptotically as: 



lim Cbc-'^Iaf{Ps) < J log 



hr 



mm 



The apparent looseness of the gaps computed above compared to the corresponding gaps in 
[|4| arises from the series of simplifications made to reduce the problem (P2) to the problem 
(P4) and the particular definition of lower bound used in Q. 

C. Type B Relay Networks 

No attenuation, constant total relay power: As claimed above, the upper-bound in Proposition]?] 
holds for Type B networks too. Therefore we only address the asymptotic gap between Cbc 



and the lower bound in (19). We have 



lim Cbc - 2/af(Ps) = I log(l + l/Q) 

The actual capacity C of the relay network in this case too is bounded by \{Cbc — ^log(l + 
l/Q)) < C < Cbc and bound gets tighter with increasing Q. 
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VI. Conclusion and Future Work 

We provide a framework to analyze the performance of the AF relay schemes in a general class 
of wireless relay networks. We demonstrate the effectiveness of the proposed framework in two 
ways. We first show that some well-established results for Gaussian relay networks (for example, 
those in Q) can be derived in a much simpler manner. Then, we show that the achievable rates 
for more complex Gaussian relay networks, which cannot be addressed with existing methods, 
can be derived straightforwardly. The proposed framework allows for the computation of tighter 
amplify-and-forward lower bounds to the capacity of Gaussian relay networks. 

Future Work: We have provided the analysis of AF scheme for different scenarios with the 
assumption of instantaneous relay operation. However, the resulting performance is not optimal 
when delay in relay operation is included, as the constant power allocation is not optimal for the 
original problem ( [TT] ). We also want to formally prove the hardness of optimization problem ( [TT] ) 



and develop efficient approximation schemes for it. An extension of our work also facilitates the 
computation of achievable rates for analog network coding scenarios for non-layered networks 
and in low to moderate SNR regimes that cannot be addressed with existing approaches. We 
plan to address it in detail in our future work. 

ACKNOWLDEGMENT 
This work has been supported by the CERG grant 412207. 



References 

[1] J. N. Laneman, D. N. C. Tse, and G. W. Womell, "Cooperative diversity in wireless networks: efficient protocols and 

outage behavior," IEEE Trans. Inform. Theory, vol. IT-50, December 2004. 
[2] Y. Zhao, R. Adve, and T. J. Lim, "Improving amplify-and-forward relay networks; optimal power allocation versus 

selection," IEEE Trans. Wireless. Comm., vol. TWC-6, August 2007. 
[3] S. Borade, L. Zheng, and R. Gallager, "Amplify-and-forward in wireless relay networks: Rate, diversity, and network size," 

IEEE Trans. Inform. Theory, vol. IT-53, October 2007. 
[4] M. Gastpar and M. Vetterli, "On the capacity of large Gaussian relay networks," IEEE Trans. Inform. Theory, vol. IT-51, 

March 2005. 

[5] K. S. Gomadam and S. A. Jafar, "Optimal relay functionality for SNR maximization in memoryless relay networks," IEEE 
JSAC, vol. 25, February 2007. 

[6] T. Cui, T. Ho, and I. Kliewer, "Memoryless relay strategies for two-way relay channels," IEEE Trans. Comm., vol. 57, 
October 2009. 



18 



[7] S. Zhang, S. C. Liew, and P. P. Lam, "Physical-Layer Network Coding," Proc. ACM MobiCom, Los Angeles, CA, September 
2006. 

[8] S. Katti, L Marie, A. Goldsmith, D. Katabi, and M. Medard, "Joint relaying and network coding in wireless networks," 

Proc. IEEE ISIT 2007, Nice, France, June 2007. 
[9] S. Katti, S. Gollakotta, and D. Katabi, "Embracing wireless interference: analog network coding," Proc. ACM SIGCOMM, 

Kyoto, Japan, August, 2007. 

[10] I. Marie, A. Goldsmith, and M. Medard, "Analog network coding in the high-SNR regime," Proc. IEEE WiNC 2010, 
Boston, MA, June 2010. 

[11] A. Argyriou and A. Pandharipande, "Cooperative protocol for analog network coding in distributed wireless networks," 

IEEE Trans. Wireless. Comm., vol. TWC-9, October 2010. 
[12] Q. You, Z. Chen, Y. Li, and B. Vucetic, "Multi-hop bi-directional relay transmission schemes using amplify-and-forward 

and analog network coding," To appear in Proc. IEEE ICC 2011, Kyoto, Japan, June 2011. 
[13] S. Zahedi, M. Mohseni, and A. El Gamal, "On the capacity of AWGN relay channels with linear relaying functions," ," 

Proc. IEEE ISIT 2004, Chicago, XL, June 2004. 
[14] R. Koetter and M. Medard, "An algebraic approach to network coding," IEEE/ACM Trans. Netw., vol. 11, October 2003. 
[15] W. Hirt and J. L. Massey, "Capacity of the discrete-time Gaussian channel with intersymbol interference," IEEE Trans. 

Inform. Theory, vol. IT-34, May 1988. 
[16] M. Avriel, W. E. Diewert, S. Schaible, and I. Zang, Generalized Concavity. Plenum Press, 1988. 

[17] S. Schaible and J. Shi, "Fractional Programming: the sum-of-ratios case," Optimization Methods and Software, vol. 18, 
April 2003. 

[18] S. Boyd, S.-J. Kim, L. Vandenberghe, and A. Hassibi, "A tutorial on Geometric Programming," Optimization and 

Engineering, vol. 8, March 2007. 
[19] M. Chiang, "Geometric Programming for communication systems," Foundations and Trends in Communications and 

Information Theory, vol. 2, August 2005. 
[20] M. Chiang, C. W. Tan, D. P. Palomar, D. O'Neill, and D. Julian, "Power control by Geometric Programming," IEEE Trans. 

Wireless. Comm., vol. TWC-6, July 2007. 
[21] M. J. Osborne, |http://www.economics.utoronto.ca/osbome/MathTutorial/QCC.HTM| 
[22] T. M. Cover and J. A. Thomas, Elements of Information Theory, lie., John Wiley & Sons, 2006. 
[23] E. Telatar, "Capacity of multi-antenna Gaussian channels," Europ. Trans. Telecommun., vol. 10, Nov./Dec. 1999. 
[24] R. G. Gallager, Information Theory and Reliable Communication. John Wiley & Sons, 1968. 

[25] B. S. Tsybakov, "Capacity of a discrete-time Gaussian channel with a filter," Probl. Peredach. Inform., vol. 6, pp. 78-82, 
1970. 

[26] W. Toms and T. Berger, "Capacity and error exponents of a channel modeled as a linear dynamic system," IEEE Trans. 

Inform. Theory, vol. IT- 19, January 1973. 
[27] L. H. Brandenburg and A. D. Wyner, "Capacity of the Guassian channel with memory: The multivariate case," Bell Syst. 

Tech. J., vol. 53, May-June 1974. 
[28] U. Grenander and G. Szego, Toeplitz Forms and Their Applications, lie,, Chelsea Publishing Co., New York, 1984. 
[29] R. M. Gray, "On the asymptotic eigenvalue distribution of Toeplitz matrices," IEEE Trans. Inform. Theory, vol. IT-18, 

November 1972. 



19 



[30] A. J. Goldsmith and M. Effros, "The capacity region of broadcast channels with intersymbol interference and colored 

Gaussian noise," IEEE Trans. Inform. Theory, vol. IT-47, January 2001. 
[31] C. Choudhuri and U. Mitra, "Capacity of relay Channels with ISI and colored Gaussian noise," Proc. IEEE ISIT 2009, 

Seoul, Korea, June 2009. 

[32] S. Agnihotri, M. Chen, and S. Jaggi, "Approximation schemes for a particular class of global optimization problems 
encountered in the capacity analysis of the ISI channels with colored noise," Under preparation, 2011. 



Appendix A 
Proof of Lemma[T] 

Following Gallager |24|, we define the information rate for the channel in ^ as 

J(P„/3) = lim lN{Ps,f3), (A.l) 

where 

/jv(P„ f3) = N-'l{x,[0, N-l];yt[0,N-l]), (A.2) 

with the assumption 

x[-D-\-l] = iO,...,0) 

It should be noted that as D'^ is finite, this assumption does not affect I(Ps,f3), but it simplifies 
the analysis as we show later. 



Following [15 1, we introduce a new channel model whose maximum achievable rate can be 
easily computed. The capacity of this channel can be, subsequently, used to compute I{Ps)- The 
output of the new channel at time instant n is given as 

N-l 

= '^hjXsiin - j) mod N] (A3) 

j=0 
N-l 

+ [(^ ~ j) ^od N] 

j=0 

N~l 

+ . . . + ^ hM,jZM[in - j) mod N] 

3=0 

+Zt[n\, <n < N, 

where N > D"^. Extending the unit-sample response of the channel for signal h[0, D^] with A^ — 
D'^ — 1 zeros, we define h[0, A^ — 1] = (ho, hi, ... , hjys, 0, 0, . . . , 0). We can similarly extend the 
corresponding unit-sample responses of the channel for noise samples zi, Z2, ■ ■ ■ , zm and define 



20 



h[0,N-l] = (/ii,o, • • • , h,,D^,0, 0, . . . , 0), h2[0,N-l] = (/i2,o, /i2,i, • • • , h,D^,0, 0, . . . , 0), 
. . ., /iAf[0, — 1] = (/iA/,0, /iA/,1, . . . , /im.da^, 0, 0, . . . , 0), respectively. With the definitions of 
these extended sample responses, we can rewrite ( |A.3[ ) symbolically as 



yt[0,N-l] = Xs[0,N -l]®h[0,N -1] (A.4) 
+2i[0,Ar- 1] ® /ii[0,Ar- 1] 
+ . . . + zm[0, N -l]®hM[0,N - 1] 
+z40,iV-l], 

where ® denotes the circular convolution operator. The transmit power constraint ([2]) leads to 
the power constraint 



E[x^i[n]] < Pi, 0<n<N, 



(A.5) 



for the new channel model. 



The information rate for the channel defined by ( |A.3[ ) and ( |A.5[ ) is defined as 

iN{Ps,f3) = N-'l{xs[0,N- l];yt[0,N- 1]) 



Taking the Discrete Fourier Transform (DFT) on both sides of (A.4) we have 



Yi = H,Xi + Hi,iZiA + ■■■ + HmaZma + Zi, 0<i<N, 



(A.6) 



(A.7) 



where Yi, Hi, Xi, Hi i, Zi j, . . ., Hm,!, Zm,!, and Zi are the components of DFT{yt[0, N — 
1]}, DFT{h[0,N -1]}, DFT{xs[0,N -1]}, DFT{h[0, N - 1]}, DFT{z,[0, N - 1]}, . . ., 
DFT{hM[0, - 1]}, and DFT{zm[0, - 1]}, respectively. 

Dividing both sides of this by Hi and transforming the resulting first [N/2\ + 1 components 



with the transformation [15. , eq. (24)], we get the following equivalent form of ( |A.3[ ) in the 
transform domain: 

r/ = x; + z;, + ... + z;,, + z:, o<z<Ar (a.8) 

where X'-, Z[^, .. ., Z'j^^^, Z[ are obtained from transform [15, eq. (24)] with Bi = X^, Bi = 
Hi^iZi^i/Hi, Bi = HM,iZM,i/Hi, and Bi = Zi/Hi, respectively. 
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Note that 



0<i< N 



where * denotes the complex conjugate. Therefore, it follows from [ 15 Lemma 1] that Z'^^, 1 < 



m < M, in ( |A.8[ ) are statistically independent zero mean Gaussian random variables with 
variance 



a. 



m,i 



Na 



\H, 



2 I -'-'r?i,i I 



0<i < N. 



(A.9) 



Similarly, Z'^ are statistically independent zero mean Gaussian random variables with variance 

1 



0<i<N 



(A.IO) 



Therefore, the equivalent transform domain channel model ( |A.8| ) is a set of parallel discrete 
memoryless additive Gaussian noise channels with total noise variance 



M 



1 + Em=l 



m=l 



0<i< N 



(A.ll) 



Further, using [15 Lemma 2], it follows that the transformed inputs < i,< N, are i.i.d. 



zero mean Gaussian random variables with variance NPs. Therefore, the equivalent transform- 



domain channel model in ( |A.8| ) is a set of parallel, independent, discrete, memoryless additive 
Gaussian noise channels with zero mean Gaussian inputs of variance NPg and total noise variance 
given in ( |A.l 1[ ). This implies that the average mutual information of i^^ component channel is 
then the capacity of point-to-point AWGN channel, given by 



/(x;,F/) = -iog 



1 + 



p. 



\Hi\ 



a' 



1 + Elli \H. 



|2 

'-m,i\ J 



(A.12) 



As the A^ component channels are mutually independent, therefore I]\f{Ps,f3) in ( |A.6| ) is given 
by 



N-l 



i=0 



1 



a' 



1 + =1 1^. 



2 

m,i I J 



(A.13) 
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Further, p?, Theorem 4] states that 

{l-D'/N)lN-Ds{Ps,fi) < /iv(Ps,/3) 

< (l + DViV)J^+^.(P„/3) 



Therefore, we finally have I{Ps,(3) in ( |A.1| ) as 

/(P„/3) = lim /;v(Ps,/3) 



(A.14) 



lim — — log 



iv^oo 2iV 

(a) ^ r 1 

27r 



j=0 



log 



1 + 



p. 



a- 



P.. 



1 + e:=i i^, 



(J ^ r 

27r 



log 



P. 



where (a) follows from the property of Riemann integrals as stated in [15 Lemma 5] and (6) 
follows from the fact that \H{\)\ = |iJ(-A)|,|A| < tt. It should be noted that H{\) is the 
transfer function of the filter with unit-sample response (ho, /ii, . . . , /i^s), given as follows 



■ijX 



(A.15) 



j=0 



Similarly, Hm{X) is the transfer function of the filter with unit-sample response {hm,o, ^m,i, ■ ■ • , hm,D^), 
1 <m < M, given by 



ijX 



-1 



(A.16) 



j=Q 

Appendix B 
Expressions for \H{X)\'^ and |iJm(A)p 



Following (|9]), we have 

\H{X)\^ = H{X)H*{X) 

D" D" 



j=0 j=0 

^ /i^ + 2( ^ hjhj+i) COS A + 2( ^ hjhj+2) cos(2A) 

j=0 j=Q j=0 

+ ... + 2hohDs cosiD'X) 



(B.l) 



+ v4i cos A + . . . + Ad" cos(D^A) 



(B.2) 
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Comparing the coefficients of cos(jA), 1 < j < -D'^, in ( |B.1[ ) and ( |B.2| ), defines Aj. 
Note that we have the following relation between AiS and hiS: 



(B.3) 



j=0 



Noting that cos(r;,A) can be expressed in terms of T„(x), the Chebyshev polynomials of the 
first kind, as cos(nA) = T„(cosA), we can rewrite \H(\)\'^ as 



\H{X)\^ = J2 (-l)'^2j + ( J2 + l)^2,+i) cosA + . . . + cos^(E.4) 

j=0 j=0 



A'q + A[cos\ + . . . + A'^s cos^ A 



(B.5) 



Comparing the coefficients of cos-' A, 1 < j < D'^, in ( |B.4[ ) and ( |B.5[ ), defines A'j. 

Similarly, following (|9]), the coefficients A'^j and ^4^,^, 1 < j < -D™, are obtained as follows: 

j=0 j=0 

= X] ^"^.i + 2( X] hm,jhm,j+i) COS A + 2( ^ hm,jhm,j+2) cos(2A) 

+ . . . + 2hmflhm,D^ C0S(D™A) 

^m,o + cos A + ... + A^,^.. cos(D™A) (B.6) 

[D™/2J [D'"/2J 

E (-l)'^2, + ( 5^ (-ir(2j + l)^,+i)cosA 
+ ... + 2^ Men. cos^ A 

^m,o + ^m,i cos A + ... + A^^^^m cos^" A (B.7) 



(a) 



The coefficients Am^i and A'^j are defined by comparing the cosine terms in (a) and (6), 
respectively, with the corresponding terms in the previous equations. 
We have the following relation between Am,iS and /im.is: 



\Hr, 



(B.8) 



1=0 
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We can rewrite J2m=i l-^ml-^)!^ 



M 



1 + 5^|//™(A)p 

71 = 1 

M 

[l + Y, ^™-o) + (Y, ^-.i) cos A + . . . + cos(D"A) 



M A/ A/ 



m=l m=l m=l 

= 5o + Si COS A + . . . + Sd« cos(D"A) (B.9) 

M M M 

= (1 + 5Z ^'m,o) + iY. ^-.l) COS A + . . . + A'^^o^) cos^" A 

771=1 m=l m=l 

= s; + s;cosA + ... + S^„cos^"A (B.IO) 

where = max(D^, . . . , D^). The coefficients i?i and are defined by comparing the cosine 
terms in (a) and (6) with the corresponding terms in the previous equations. 

Appendix C 

Illustration of the Problem P2 for Two Types of Relay Networks 

Type A Relay Network: Let us first consider the relay network configuration in Figure [2} This 
is the simplest instance of the general class of relay networks we consider in this paper. For 
M = 2 relay nodes, the Type A relay network is given in the Figure |7| We illustrate the problem 
formulation ( [TT] ) for the network in Figure |7] 
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The ISI channel between the source and the destination for this network is: 

y[n] = hstx[n] + {hsif3ihu + hs2(^2h2t)x[n - 1] 
+(5ihitZi[n - 1] + (i2h2tZ2[n - 1] + z[n], 



(C.l) 



where 



< < al 



Pi 



< < 



P2 



2,max i^2^p^ _^ 



Comparing ( |C.1| ) with Q defines corresponding modified channel parameters as follows. 

ho = hst, hi = hsiPihit + hs2(32h2t 
hi.i = (3ihu, h2A = (32h2t 



Using definitions in ( |B.2[ ) and ( |B.9[ ) results in: 

Ao = hl + hi, Ai = 2hohi 
Bo = l + hl^ + hj^ 



Therefore, we have the following formulation of the problem ( [TT] ) for the relay network in 
Figure |7| 



Iaf{Ps) = max 



2n 



log 



1 + 



subject to 



Ps Ao + Ai cos A 
^ Bo 

< Pi 

< 



Pt{hi,P, + a') < Pi 
f^M^Ps + a') < P2 



dX 



(C.2) 



Type B Relay Network: Let us next consider the network in Figure [3] For M = 2 relay nodes, 
the Type B relay network is given in the Figure [8] For M = 2 relay nodes, the source-destination 
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channel for this relay configuration is: 

y[n] = hstx[n] + +z[n] 

+ {hsif3ihu + hs2l^2h2t)x[n - 1] 
+ {hsi(3ihi2(32h2t + hs2(32h2i(3ihit)x[n - 2] 
+f3ihuzi[n - 1] + f3ihi2f32h2tzi[n - 2] 
+f32h2tZ2[n - 1] + f32h2il3ihitZ2[n - 2], 



(C.3) 



where 



Pi 



o<Pi< a 



{hi,mi + hi,)Ps + {mi + 1)^' 

P2 



2,max 



Comparing ( |C.1[ ) with (|7]) defines corresponding modified channel parameters as follows. 

ho = hst, hi = hsiPihit + hs2(32h2t, h2 = /i^ i A 12/32/^2* + h,2P2h2iPihu 
hi,i = Pihit, ^1,2 = (3ihi2(32h2t, /i2,i = (32h2t, ^^2,2 = P2h2iPihit 



Using definitions in ( |B.2[ ) and ( |B.9[ ) results in: 

Aq = hl + hj + hi, Ai = 2{ho + h2)h, A2 = 2/io^2 

Bo = 1 + /l^ 1 + 2 + ^2 1 + ^2 2 > ^1 = 2/il,l/il,2 + 2/^2,1^2,2 
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Therefore, we have the following formulation of the problem ( [TT| ) for the relay network in 
Figure [8| 



Iaf{Ps 



max 



27r 



log 



1 + 



subject to 



Ps Y.i=Q A cosiA " 

(T^ Bq + i?i COS A 

< Pi 
< Pi 



dX 



Pl{&lhl, + hl,)P, + {Plhl, + l)a') < P, 
Note the nonconvex nature of the constraint set. 



(C.4) 



Appendix D 

Proofs of Lemmas and Propositions in SectionHV] 
Proof of Proposition^ We have from (|5]) the following nonzero modified channel parameters: 



M 



hm,i = /9, for all m, 1 < m < M 



Therefore, from (B.3), (B.6), and (B.9), we have: 



Ao + Ai = H + Mfif 
Bo = l + M/32 

Substituting these in the problem (P4) establishes the Proposition. 

Proof of Lemma |2]- Setting the first-derivative of the objective function equal to zero gives 
/3 = 1 as globally optimal solution. 

Also, the objective function is quasiconcave over the domain of optimization. Note that the 
objective function is nondecreasing for all /3 < 1 and nonincreasing for all /3 > 1, then its 
quasiconcavity follows from pTj Proposition 3]. 

Proof of Proposition^ We have from (|5]) the following nonzero modified channel parameters: 



ho = hst, hi = (3 hsihj 



i=l 



hm 1 = Phmt, for all m, 1 < m < M 
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Therefore, from (IB. 31), (IB. 61), and (IB. 91), we have: 



A/ 



M 



i=l 



Substituting these in the problem (P4) establishes the Proposition. 

Proof of Lemma ^ Setting the first-derivative of the objective function equal to zero gives 
Popt = Z^'^m" .2 as the globally optimal solution. 

Noting that the objective function is nondecreasing for all (3 < (3opt and nonincreasing for all 



(3 > f3opt, establishes its quasiconcavity from [21 Proposition 3] 



Proof of Proposition^ We have from (|5]) the following nonzero modified channel parameters: 

M-{i~l) 

hi= J2 = {M -{i- l<i<M, 

hi,i = f3\l <i < M, 
h2,^ = (3\l<i<M -1, 



Therefore, from ( |B.3| ) and (|B.8|), we have: 



\Hrr 



(3-1 (3-1 ^' {(3 - 1)2 
1 < m < M 



1 + /3M 



2, pM-{i^l)_^^ 2 



Substituting these expressions in (14) establishes the proposition. 



Appendix E 
Proof of Proposition |4] 

Broadcast upper bound: For the broadcast cut, we have the following maximization problem: 

max I{Xs] Yt,Yi, . . . ,Ym\Xi, . . . , Xm) 
subject to the power constraints 

EX^ < Ps and EXf < i^, 1 < z < M 
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As the destination node knows the channel gains hit,i G St, we have 

/(X,; Yt, Fi, . . . , Ym\Xu ...,Xm) = /(X,; Yt, YgJ 

where Yt = Yt - T^^^St ^n^i^ = {Yi = Yi - Y.j&M{i)^sh3i^j ■ ^ ^ ^s}- Therefore, the 
maximum achievable rate of information transfer across the particular broadcast cut Cbc of the 
relay network is given by 

Cbc= max I{Xs;Yt,YsJ 

PXs-EXl<Ps 

The right-hand side can be evaluated using the capacity formula for Gaussian vector channels with 



fixed transfer function in [23 1. We have a channel with scalar input and vector output, therefore 
corresponding channel matrix has only one singular value given by a(Ns) = h'j.^ + J2i£Ss ^"si- 
Hence, we have 



Cbc = log 



1 + ^a{Ns) 



(E.l) 



Multiple-access upper bound: For the multiple-access cut, we have the following maximiza- 
tion problem: 

max I{Xs, Xi,..., Xm; Yt) 

subject to the power constraints 

EX^ < Ps and EXf <Pi,l<i<M 
However, given that Yt = hstX^ + T^ieXs^ ^a^i^ we have 

J(X„ Xi, . . . , Xm; Yt) = /(X„ Yt) 
Therefore the maximization problem above reduces to 

max I{Xs,Xs,;Yt) (E.l) 

subject to the power constraints 

EX^ < Ps and EXf < PiJ e St (E.3) 
Let us relax this power constraint to 

EX^ + J2 EXf <Ps + J2Pi (E.4) 

i&St i&St 
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The solution of the problem defined by ( |E.2[ ) and ( |E.3| ) cannot be larger than the solution of the 
problem defined by ( |E.2[ ) and ( |E.4| ). The latter can be evaluated by using the abovementioned 
result in [ |23| again. We have a channel with vector input and scalar output, therefore the 
corresponding channel matrix has only one singular value given by d(Nt) = h^^ + J2ieSt ^'it- 
Therefore, the maximum achievable rate of information transfer across the particular multiple- 
access cut Cm AC of the relay network is given by 



C 



MAC 



log 



1 + 



-d{Nt 



(E.5) 



As the capacity C of the relay network must be smaller than the minimum of Cbc and Cmac^ 
this proves the proposition. 



